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A general method, an algorithm, and results of computer modeling of destruction of hardened sheet glass un¬ 
der conditions approximating testing of real products are described. Recommendations for further improve¬ 
ment of the program are given. 


Computational methods [1,2] for predicting the nature 
of destruction of hardened sheet glass make it possible to de¬ 
termine the total number of fragments N per prescribed sur¬ 
face area of the article (5x5 cm), assuming uniform distri¬ 
bution of them on the surface and identical shape 
(parallelepiped) and size of each fragment. Testing of sam¬ 
ples according to GOST 5727-88 yields a certain spread in 
these parameters. 

The present paper attempts to give a computerized 
graphical representation of destruction of sheet glass under 
conditions approximating real conditions for product testing. 

The development is based on analytical dependences ob¬ 
tained earlier [2]. 

In particular, the average length of a side of the base of a 
fragment is equal to [2] 
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where E is the modulus of elongation; y is the specific sur¬ 
face energy; ct is the central tensile stress; v is the Poisson 
coefficient; 5 is the glass thickness. 

The average number of fragments for the length B of a 
side of a square plate will be 

" = B 2 /Zl- 

It is assumed that the crack grid is formed by N r rows, 
each containing N r fragments. The destruction representation 
scale is determined: 


mas = 400/5, 

where 400 is the sheet side size in pixels (a pixel is a conven¬ 
tional-measurement unit for a computer display; one conven¬ 
tional unit ranges from 0.25 to 0.28 mm). 


The average length of a base side of a fragment repre¬ 
sented on the computer display is 


Z av = (B/N r ) mas. 


Experimental data show that fragments become larger as 
they recede from the center of destruction (by the standard, 
this is the geometric center of the article). We assume that 
within the confidence interval ± 20%, the variation in Z (Z is 
the length of a base side of a fragment represented on the 
computer display) as it recedes from the center of destruction 
will be 


0.8Z av < Z < 


.2Z„ 


These initial data were used as the basis for an original 
computer program “SETGRAPH” providing a graphical rep¬ 
resentation of destruction of hardened glass. The general 
structure of the program is as follows. 

One of the main purposes of the program is to position 
the center of each fragment, taking into account the preset 
variation in the length of a fragment side. First, the centers of 
fragments of the first row are positioned, the ordinal number 
of the fragment k ranging from 1 to N r . The distance between 
fragment centers is found from the formula 


AAD [wi, k ] 



B, + 0.8Z.„ 


where m is the ordinal number of the row in which the frag¬ 
ment is positioned on the display; k is the ordinal number of 
the fragment in the row; 0.8Z av is the arbitrarily adopted mi¬ 
nimum length of a fragment side; B ( is the absolute value of 
the difference in the length of a side of the preceding and fol¬ 
lowing fragments in the row with the number /«, 
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where 0.4 Z is the difference between the maximum and 

<1V 

minimum length of a fragment side. 

As a result, the coordinate along the X axis for the center 
of fragment A in row m will be 

X0 [m, k] = X 0 [m, k - 1 ] + AAO [m, k ]. 

The value of the coordinate along the Y axis for the cen¬ 
ters of the fragments in this row is constant. 

Next we construct the angles between the base sides of 
each fragment shown on the display (Fig. 1). A fragment is 
arbitrarily split into four sections, the coordinates of each 
vertex being determined: 

X (m, k , i ) and Y (m, A, i ), 

where i is the number of the fragment angle, ranging from 1 
to 4. 

Let us assume that all vertices of the fragment base are 
equally far from its coordinate origin (Fig. 1), i.e., are posi¬ 
tioned on a circumferences with the radius /?[/«, A], The 
value of the radius is different for each fragment and is found 
as follows: 



Fig. 1. Scheme of the glass fragment structure: /- IV) quarters into 
which the fragment is split. 
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R [/«, A] = Z [in, A'] 
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where 


Z A ] = 




A + m 
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We assume that the shortest distance from a base vertex 
to the Y axis is 

L x [m, A, / ] = random (0.6/?) + 0.335/?, 



where 0.335 R is the arbitrarily adopted minimum possible 
distance from the vertex to the Y axis; 0.6 R is the difference 
between the arbitrarily adopted maximum (0.935 R ) and mi¬ 
nimum distances from the vertex to the Y axis. 

The symbol “random(0.6/?)” in the above relationship is 
a function that generates a random number in the interval 
( 0 . 6 /?). 

The shortest distance from a base vertex to the X axis is 
L y [m, A. i ] = ^Jr 2 [iii, k]- Lr x [m, A, /]. 


Fig. 2. Example of the graphical representation of destruction of 
hardened glass (AA = 81, size of a sample side 50 mm). 

ing from 0 to 200. An example of the graphical representa¬ 
tion of destruction of hardened glass is shown in Fig. 2. 

A special feature of the developed program is the possi¬ 
bility of introducing any distribution law for the size of a 
fragment side, for instance, the Gauss law, and of changinng 
the shape of the simulated fragments. 


The base vertices on the display are connected by conti¬ 
nuous lines, which represent the crack grid of the destroyed 
hardened glass. 

The initial data (processed by the SETGRAPH program) 
make it possible to obtain a crack grid for hardened sheet 
glass up to 0.9 m 2 in size with the number of fragments rang- 
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